7r-phases in balanced fermionic superfluids in spin-dependent optical lattices 
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We study a balanced two-component system of ultracold fermions in one dimension with attractive 
interactions and subject to a spin-dependent optical lattice potential of opposite sign for the two 
components. We find states with different types of modulated pairing order parameters which 
are conceptually similar to 7r-phases discussed for superconductor-ferromagnet heterostructures. 
Increasing the lattice depth induces sharp transitions between states of different parity. While the 
origin of the order paramter oscillations is similar to the FFLO phase for paired states with spin 
imbalance, the current system is intrinsically stable to phase separation. We discuss experimental 
requirements for creating and probing these novel phases. 
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One of the most intriguing examples of the interplay 
of superconductivity and magnetism is the Fulde-Ferrel- 
Larkin-Ovchinnikov (FFLO) phase, where Zeeman split- 
ting of the Fermi surfaces is expected to lead to spatial 
oscillations of the pairing amplitude. It is difficult to 
obtain such phases in superconductors, since the orbital 
effect of the magnetic field is typically much larger than 
the spin Zeeman splitting. Several proposals have been 
made, however they remain controversial pQ. For exam- 
ple, FFLO phase has been discussed in the context of 
heavy fermion CeCoIn5 superconductors [3], but al- 
ternative interpretation in terms of competing magnetic 
order has also been given [4|. So far the only unam- 
biguous demonstration of FFLO-like physics has been 
achieved in heterostructures of ferromagnetic and super- 
conducting (F/SC) layers [5], where proximity coupling 
through ferromagnetic layers results in superconducting 
7r-junctions (see [1] for a review). We note that 7r-phases 
arising from a different mechanism than FFLO have also 
been discussed for high-T c cuprates [TJ [5] . 

Recently, in cold atoms, there has been a large body 
of work, both experimental and theoretical, aimed at 
achieving FFLO states. The biggest difficulty is that 
FFLO phases are fragile and extremely susceptible to 
phase separation and the experimental situation remains 
unclear [71112). In this paper we propose a novel system 
of ultracold fermions in an optical lattice O [H] which 
can be used to observe FFLO type states with oscillat- 
ing pairing amplitude. The system which we discuss is 
somewhat similar to F/SC heterostructures and should 
be stable against phase separation. Our proposal relies 
on the ability to create spin dependent optical lattices 
and we find that beyond a certain critical strength of 
such optical potential, the superconducting pairing am- 
plitude becomes a sign changing function (we will refer 
to such states as as 7r-phases, see Fig. [2]). 

The gap profile in the ground-state depends on the 
wavelength of the lattice, A, the strength of the poten- 
tial, Vq, and the interaction strength. In Fig. [T] we present 




FIG. 1: Phase diagram showing the emergence of 7r-phases 
for a spin-dependent lattice potential of wavelength A and 
strength Vo for interaction strength giokF / Ef ~ —2.04 
and zero temperature. A gradient of colors gives A max := 
max|Amj, the largest Fourier component amplitude of the 
gap. The black lines indicate transitions from gap profiles 
with zero (BCS), two (n), four (2tt) and six (3tt) zero-crossings 
per unit cell. The dashed FF line and the BCS-FF arrow are 
from a Fulde-Ferrell calculation in the homogeneuos system 
and is explained in the text. 

the (Vo,A) phase diagram showing the transitions from 
constant gap to the 7r-phases with several zero-crossings 
per unit cell in the pairing amplitude. A color gradient 
gives the largest Fourier component amplitude of the gap 
(see below) and the black lines indicate the transitions. 
We clearly see 7r-phases occurring in a broad range of 
A restricted from below only by the coherence length as 
we will discuss. The emergence of oscillations in the gap 
gives clear signatures in the Fourier transform. We will 
demonstrate how the rapid-ramp techniques can be used 
to observe these states in time-of-ffight measurements. 
We also suggest ways to make spin dependent large wave- 
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length lattice potentials in the high-field regime as is 
needed to access 7r-phases. 

The quasi-lD system we study is described by the ef- 
fective Hamiltonian [TJ] 

H - - = 
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where gip, is the effective ID coupling constant. We use 
g\p>kp/Ep = —2.04 as in [T5] corresponding to an inter- 
action to kinetic energy ratio of —rngm/fi 2 = 1-6. This 
is neither weak coupling nor the unitary limit. 

We consider a balanced system but introduce chem- 
ical potentials fi a since the optical lattice potential is 
spin-dependent. For the main part of this work we use 
Vf(x) — Vq cos(27rx/A) and Vj,(x) = —Vf(x) in which 
case fj^ = jttj,. In the non-interacting system, the spin- 
dependent lattice spatially displaces the degenerate so- 
lution of the two components as V) is increased. In a 
simple-minded picture, pairing of these states will gener- 
ate spatial variation in the order parameter which is the 
origin of 7r-phases. 

We solve the Hamiltonian of Eq. ([!]) in mean- 
field theory by using the inhomogeneous Bogoliubov- 
deGennes (BdG) ansatz for the field operator 4 r (7 (x, t) = 
Y.kl u k'y( x ) e ~ lu ""' tc kcr + crvks(x)e lu "" >t cl (f }, where the c 
and c' denote the quasiparticles and the sum runs over 
> with k the quasiparticle index composed of a 
quasimomentum and the band index. The mean-field 
equations are 
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where n a (x) = (^(x)^ a (x)) and A(x) = 
~9id(^i(x)^i(x)) . These equations can be solved 
self-consistently for densities and gap through 
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where f(ujka) = 1/(1 + exp(hujka/kBT)) is the Fermi- 
Dirac distribution. In Eq. ^ we use the periodicity of 



the optical lattice to do a Fourier decomposition. Notice 
that these equations only contain it, v for a =t and the 
sums over uj k ^ are unrestricted [15]. We have explicitly 
checked the convergence of our numerical solutions by 
extending the cut-off on the basis size. 

The mean-field BdG ansatz does not take into account 
soft collective modes of the order parameter which, in 
principle, lead to the power law decay of the supercon- 
ducting correlations. However, the BdG approach de- 
scribes the ground state energy in our parameter regime 
well, which is determined by correlations on the scale of 
the BCS correlation length, £ = hv F /A [13 US]. Thus 
we expect it to also correctly capture the competition 
between 0- and 7r-phases. 

In Fig. [2]we show (n-[(x) — ni(x))/kp and A(x)/Ep as 
functions of 2x/X with fc^A = 30 for amplitudes Vo/Ep = 
0.39, 0.40, 0.77 and 0.78. Here we notice a sudden change 
in the gap profile from an even to an odd function (around 
x = 0) at the definite value Vo/Ep = 0.39, and again at 
Vo/Ep — 0.78. We will give arguments as to why this 
occurs in the following sections. The signature of the 
new phases are even more clear in Fig. [3] which shows 
the largest components of |Am|. Here we see a very clear 
jump between even and odd Fourier components as Vq 
is increased. This transition constitutes the main result 
of our paper and below we propose a way to observe 
the 7r-phases which is clearly distinguishable from other 
oscillatory behaviors in the gap. 

The spin-dependent lattice we use here is invariant un- 
der spatial reflection. In addition, for the balanced sys- 
tem we study, the symmetry of the lattice potential im- 
plies that the densities are even and interchanged every 
half wavelength (A/2). We can use this observation to 
restrict the functional form of the gap. In the absence 
of any currents, the gap obeys A(x + A/2) = ±A(x). 
Combined with the full periodicity A(x + A) = A(x), we 
see that either only even or only odd Fourier components 
survive which facilitates its unmistakable detection. 

Our spin-dependent lattice potential effectively acts as 
a spatially varying magnetic field. In contrast, in homo- 
geneous systems the wavevector of a Fulde-Ferrell (FF) 
state increases monotonically with increasing chemical 
potential difference. Using the results presented in [15] . 
we determine which chemical potential difference Vq (at 
fixed density) is needed so that the order parameter of 
the FF state has the wavevector 2-7r/A corresponding to 
the lattice used in our calculation. The result is pre- 
sented as a dashed line in Fig. [T] This line indicates 
where FF solutions exists, however only for Vq > 0.27 
(the BCS-FF arrow in Fig. [IJ is it a global energy min- 
imum. The qualitative agreement clearly demonstrates 
the connection to FFLO physics. The difference for our 
system is that the wavelength of variations in the gap 
does not vary continuously but changes at discrete val- 
ues of Vq since the oscillations must be commensurate 
with the lattice potential. FFLO states are similar to 
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FIG. 2: The polarization (rif(x) — ni{x))/k F in dot-dashed blue and pairing A(x)/E F in thick red as functions of 2x/X, for 
the case of 0-phase (left), 7r-phase (two figures in the middle), and 2n- phase (right) at zero temperature. The dashed and 
dotted black lines in the left plot show the spin-dependent lattice potential. Here k F X — 30, T = 0, giok F / E F = —2.04, and 
V t (x) = -Vi(x) = V cos(27rxJX). 
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FIG. 3: Plot of the absolute value of the Fourier components 
of A(a;) in Eq. (J3j) as function of the lattice potential strength 
Vo (parameters as in Fig. |2j. Only even components are non- 
zero for Vo/E F < 0.38, whereas only odd ones are non-zero 
for 0.39 < V /E F < 0.77, and so forth. This is the tell-tale 
sign of the transition from the 0- to the 7r-phase. 



7r-phases since modulation of the pairing field generates 
a lower energy solution. However, the present proposal 
differs from FFLO since we do not have a global spin 
imbalance and the system is intrinsically stable to phase 
separation. 

The transition between the different 7r-phases can be 
explained by an energy balance argument. They are 
driven by the competition between interaction (pairing 
and Hartree terms) and potential energy in the lattice. 
First consider the situation where the gap and densities 
are almost constant and even functions, A e (x) = Ao 
and n t (x) = n^x) = n /2 (V Q /E F < 0.38 in Fig [2), 
and contrast this with the situation where the gap is 
an odd function. Let us for simplicity assume the same 
magnitude of the gap and introduce a corresponding os- 
cillation in the densities, thus A°(x) = Ao sin(27rx/A) 
and n°^(x) = n /2 =p (Sn/2) cos(2nx/X). In the long- 
wavelength limit we can neglect the kinetic energy and 
the energy densities of the even and odd state can be 



= dx 
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where L is the system size. From our ansatz we get 
v e = gi D nl/4:+Al/g 1D and v Q = g XD nl/A+Al/{2g 1D )- 
Vo5n/2 — (gi F) /2)(Sn/2) 2 . If we determine the density 
variation in the odd state by requiring minimal energy, 
we find v = gi D nl/A + &l/(2g 1D ) + V 2 /(2g 1D ) from 
which it follows that the constant even solution is lower in 
energy until V = A . Taking A /E F - 0.3 from Fig. [2] 
gives V /E F ~ 0.3. Numerically we find V /E F - 0.39, 
about 30% higher. At small A we expect large deviations 
from this estimate. This is caused by the neglected ki- 
netic term that grows with decreasing A and pushes the 
jump to larger Vq. Non-commensurate solutions should 
be thermodynamically unstable by the same argument as 
they do not benefit from lowering of the potential energy. 
Initially we did search for non-commensurate solutions 
but as expected we found only commensurate ones. 

The transitions we find are very sharp as illustrated in 
Figs. [2] and [3] The 0- and 7r-phase have different parities 
and we therefore have a crossing of ground-states as we 
tune Vq. We test the stability of our predictions by us- 
ing a potential that has Vf(x) = — 3Vj,(a;)/2. As Fig. [4] 
shows, a sharp transition occurs also in this case. Even 
though this potential breaks A/2 symmetry, there is still 
conservation of parity, thus A_m = iA^, and a sharp 
transition still occurs. However, now the order parameter 
constains both even and odd Fourier components. 

The results presented in Figs. [3] and [4] have fc^A = 30. 
For smaller A ~ £, the lattice drives the system into the 
normal state before showing any noticeable oscillations 
of the gap. For the interaction strength gmk F / E F — 
—2.04, we estimate that k F X > 12 is necessary to support 
observable 7r-phases. In the phase diagram in Fig. [I] the 
suppression of the gap at the transitions at small A is 
clearly seen. Fig. [T] also demonstrates that more zeros of 
the gap per unit cell could be accessible in experiments. 
The presented results are for the zero temperature case. 
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FIG. 4: Same as Fig. [5] but with = -3Vj.(a;)/2 = 

Vb cos(27ra;/A). 

For the parameter regime in Fig. [I] we find that the phase 
diagram is qualitatively the same up to T ~ O.lTp. We 
have not performed a detailed analysis of the transition 
to the normal state outside of this parameter regime or 
for even higher temperatures. An interesting aspect of 
this system is that the spin-dependent potential leads to 
an appreciable triplet component in the pairing [T7] . 

In order to detect the 7r-phase, the rapid ramp tech- 
nique can be used to transfer opposite spin pairs into 
molecules on the Bose-Einstein condensate (BEC) side 
of the Feshbach resonance as shown in [T8Tf2"0] . We 
can calculate the number of molecules with center-of- 
mass momentum q := 27rm/A, with rh an integer, by 
projection of the state before the ramp into the bound 
s-wave molecular state f s {r) ~ y/2/a s exp(— r/a s )/r, 
where r is the relative distance and a s the scattering 
length on the BEC side (2TJ [22]. Assuming a s ,a± <C £, 
with aj_ the transverse confinement length, we obtain 
rig ~ (8La 2 _/g 2 D a s )K[a±/a s ] 2 \Am\ 2 Additional 
density terms give only a featureless continuous contri- 
bution to rig. For V a (x) = — V-„(x), 7r-phases are linked 
to either only odd or only even Fourier components of 
the order parameter. After the rapid ramp process, the 
BCS to 7r-phase transition is seen as the disappearance 
of a single spot and the apperance of two new spots, 
ruling out other types of oscillation in A. In contrast, 
phase-separated densities located in the minima of the 
potentials would require many Fourier components. Ob- 
servation of sharp peaks in the distribution is therefore 
a clear sign of 7r-phases. 

For experimental realization we focus here on 40 K [21] . 
6 Li is another possible candidate, although we note that 
a spin-dependent lattice is harder to implement [25] . We 
assume a ID geometry of tubes that are optically trapped 
with a superposed magnetic field to control the interac- 
tion via the Feshbach resonance at B — 202. 1G. Using 
N ~ 100 per tube of length L ~ 40/^m, we have n ~ 2.5 
parr 1 and kp = nn/2 ~ 3.93 /^m _1 . For simplicity we 



neglect the external confinement. With a± = 60.3nm 
and a s = 132. 3nm from [24], and using A.±/Ep = 0.15 
and giokp/Ep — —2.04, we find n 27r /\ — 1.58. How- 
ever, the number of tubes in [23] was 4900 so we expect 
a signal of 7742 molecules, which should be detectable. 

In the unpolarized system, the appearance of 7r-phases 
in the order parameter requires a spin-dependent lattice 
potential with a wavelength longer than the coherence 
length; A > £. To fulfil both requirements multiple lasers 
should be used [2SH2H]- To get spin-dependence there 
are several proposals and we focus on the one of [2"5] , 
The splitting is controlled by the difference in laser inten- 
sity and phase of left-circular and right-circular polarized 
light. Furthermore the transition is between the 2 S\/2 
and 2 Pi/2j 2 P3/2 imes with optical wavelength X opt ~ 770 
nm. To change A in the lattice one changes the angle be- 
tween laser and tubes. The magnetic field is aligned par- 
allel to the lasers. If 9 is the angle between the tubes and 
these lasers, the lattice wave-length is A = A opt /2cos(6 l ). 
kpX = 30 translates to 7.6 /im and 6 ~ 87°. Since this 
is almost perpendicular to the ID tube, the heating will 
also be reduced as most of the recoil is absorbed in the 
confining potential. 

In an actual experiment it has been suggested that ID 
tubes in an intermediate strength 2D optical lattice will 
give the best conditions for observing the FFLO state in 
ID and we expect this to hold for our 7r-phases as 
well. Our proposal differs from other studies on FFLO 
states in cold atom since we use an unpolarized gas. The 
Fermi surfaces are therefore identical for the two spins 
and the non-trivial pairing properties of the system are 
entirely due to the spin-dependent lattice potential. 

We thank R. Scnsarma, D. Pekker, L. Fritz, D. Weld, 
M. A. Cazalilla and F. Sols for helpful discussions. The 
authors acknowledge support from Real Colegio Com- 
plutense en Harvard, MEC (Spain) grant FIS2007-65723, 
the German Research Foundation grant WU 609/1-1, 
the Villum Kann Rasmussen foundation, CUA, DARPA, 
MURI, and NSF grant DMR-0705472. 



[1] A. I. Buzdin, Rev. Mod. Phys. 77, 935 (2005). 

[2] H. A. Radovan et al, Nature 425, 51 (2003). 

[3] A. Bianchi, R. Movshovich, C. Capan, P. G. Pagliuso, 

and J. L. Sarrao, Phys. Rev. Lett. 91, 187004 (2003). 
[4] F. Ronning et al, Phys. Rev. B 71, 104528 (2005). 
[5] P. H. Barsic, O. T. Vails and K. Halterman, Phys. Rev. 

B, 75, 104502 (2007). 
[6] E. Berg et al, New J. Phys. 11, 115004 (2009). 
[7] M.W. Zwierlein et al., Science 311, 492 (2006); Nature 

(London) 442, 54 (2006) 
[8] Y. Shin, M. W. Zwierlein, C. H. Schunck, A. Schirotzek, 

and W. Ketterle, Phys. Rev. Lett. 97, 030401 (2006). 
[9] C.H. Schunck et al., Science 316, 867 (2007). 
[10] G.B. Partridge et al, Science 311, 503 (2006); Phys. Rev. 

Lett. 97, 190407 (2006). 



5 



[11] P.F. Bedaque, H. Caldas, and G. Rupak, Phys. Rev. Lett. 

91, 247002 (2003). 

[12] Y. Liao et al, |arXiv :0912.0092v2. 

[13] R. Jordens et al, Nature (London) 455, 204 (2008). 
[14] U. Schneider et al, Science 322, 1520 (2008). 
[15] X.-J. Liu, H. Hu and P. D. Drummond, Phys. Rev. A 76, 
043605 (2007). 

[16] H. Hu, X.-J. Liu, and P. D. Drummond, Phys. Rev. Lett. 

98, 070403 (2007). 
[17] I. Zapata et al., in preparation. 

[18] M. Greiner, C.A. Regal, and D.S. Jin, Nature (London) 

426, 537 (2003). 
[19] C.A. Regal, M. Greiner, and D.S. Jin, Phys. Rev. Lett. 

92, 040403 (2004). 

[20] M.W. Zwierlein, C.A. Stan, C.H. Schunck, S.M.F. Rau- 
pach, A.J. Kerman, and W. Ketterle, Phys. Rev. Lett. 
92, 120403 (2004). 



[21] R.B. Diener and T.L. Ho, |arXiv: cond-mat /04045 17| 
[22] E. Altman and A. Vishwanath, Phys. Rev. Lett. 95, 
110404 (2005). 

[23] K [a] = r dx [°° d y^-Vf+^y'-y 2 /*) 

[24] H. Moritz, T. Stoferle, K. Giienter, M. Kohl, and T. 

Esslinger, Phys. Rev. Lett. 94, 210401 (2005). 
[25] J.K. Chin et al., Nature (London) 443, 961 (2006). 
[26] D. Jaksch, H. J. Briegel, J. I. Cirac, C. W. Gardiner, and 

P. Zoller, Phys. Rev. Lett. 82, 1975 (1999). 
[27] G. Grynberg and C. Robilliard, Phys. Rep. 355, 335 

(2001). 

[28] W.V. Liu, F. Wilczek, and P. Zoller, Phys. Rev. A 70, 
033603 (2004). 

[29] M.M. Parish, S.K. Baur, E.J. Mueller, and D.A. Huse, 
Phys. Rev. Lett. 99, 250403 (2007). 



